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We present an exploratory study that suggests that Efimov physics, a leading research 
theme in few-body quantum physics, can also induce stable many-body ground states whose 
building blocks are universal clusters. We identify a range of parameters in a mass-and- 
density imbalanced two-species fermionic mixture for which the ground state is a gas of 
Efimov-related universal trimers. An explicit calculation of the trimer-trimer interaction 
reveals that the trimer phase is an SU(3) Fermi liquid stable against recombination losses. 
We propose to experimentally observe this phase in a fermionic 6 Li- 53 Cr mixture. 


* Author contributions: SE carried out the Born-Oppenhaimer calculation and the analysis of the experimental 
feasibility of the theoretical model. AMG initiated the project and wrote a first draft of the manuscript. PN 


computed the effective trimer-trimer potential and the scattering length by the resonating group method. 
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I. INTRODUCTION 


The Efimov effect [it3], predicted by Vitaly Efimov more than forty years ago, states that 
three particles interacting through a sufficiently strong two-body interaction can occupy a tower 
of three-body bound states usually referred to as Efimov trimers. Intriguingly, these states can 
exist even if there are no two-body bound states, and feature discrete scale invariance. The 
binding mechanism is an effective three-body attraction that has no classical analogue. The 
landmark experimental observation of an Efimov trimer in a gas of Cesium atoms [4] provided 
new impetus to this problem. Theoretical features such as the discrete geometric scaling of the 
energy levels corresponding to trimers [5-8] and the universality of Efimov physics [3, El were 
observed experimentally in the vicinity of Feshbach resonances. All these experiments were 
carried out in the dilute limit where many-body effects are negligible. 

Some effects of a many-body environment on Efimov physics have already been studied 


theoretically 1 OH 171]. Pauli 


blocking due to a static Fermi sea tends to suppress the Efimov 


trimers as density increases [10—13]. By contrast, the presence of a single species Fermi sea 


enhances the stability of a single p-wave universal (non-Efimov) trimer composed of two species 
of fermions [14]. However, the more fundamental question of whether Efimov physics can 
induce novel many-body phases is only starting to be addressed in the literature. The calculation 


of the third virial coefficient in a mass-imbalanced two-component Fermi gas H and Bose gas 
at unitarity [Q suggests an Efimov-driven phase transition. For a three-component Fermi gas 
the possible existence of a trimer phase was proposed in Ref. 017 1. However, in ultracold atom 
systems, the formation of Efimov trimers entails three-body recombination processes leading 
to a significant atom number loss [3-jj], and it is unclear whether it is possible to observe this 
phase in these systems. 

Here, we propose a system closely related to Efimov physics, where a many-body trimer 
phase can be stably realized. We show that in a mass-and-density imbalanced two-component 
Fermi gas, a strong enough interaction between the two components leads to a many-body phase 
composed of trimers. These trimers are not Efimov trimers but universal trimers, also known 
as the Kartavtsev-Malykh trimers [18], which are precursors of the Efimov effect occurring at 
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large mass ratios. We explicitly investigate the interaction between the trimers to characterize 
the nature of the trimer phase. We find the trimer-timer interaction is repulsive at low energy, 
which suggests that the trimer phase is a stable Fermi liquid. 


II. TRIMERS IN MASS AND DENSITY IMBALANCED FERMI MIXTURES 


We study a homogeneous two-species Fermi mixture where the mass and number density 
of each component is denoted by and respectively. The interaction between different 
components is assumed to be attractive and near-resonant, i.e. it is characterised by a large 
s-wavc scattering length |a| » r 0 , where r 0 is the typical range of the interaction. It can thus 
be modeled by a zero-range contact interaction with scattering length a. The inverse scattering 
length 1/a goes from — oo to +oo as the strength of the interaction increases, and is zero exactly 
at the resonance. 

For equal masses rrif = m 4 = m, and equal densities rif = /q = n, the ground state evolves 
smoothly from a Bardeen-Cooper-Schrieffer (BCS) phase characterised by overlapping Cooper 
pairs in the weak-coupling limit n 1//3 a — » CT, into a Bose-Einstein Condensate (BEC) of dimers 


with a binding energy E^ = — K 2 /ma 2 in t 


le strong-coupling limit n 1//3 a —» 0 + 


1911. For 


n \ 77,4. the phase diagram is richer [|20l. 21], but the ground state is still a superposition of 


2211230, trimers and tetramers can be 


Cooper pairs or dimers rather than of trimers or tetramers. 

By contrast, for a mass-imbalanced system [2, 

formed due to the Efimov attraction developing in odd angular momentum channels. For a three- 
body system of ]”|\|, particles, Efimov trimers exist for a sufficiently large mass ratio m^/m^ > 


13.6 [2, 


23]. These trimers cannot be realized stab 


in cold-atom experiments due to their 


a : 


at short distances of the order of r 0 , 


internal instability against three-body recombination 13- 
two atoms of the trimer can recombine into a tightly bound diatomic molecule, whose binding 
energy is released by the third atom. On the other hand, for a mass ratio 8.2 < m^/m± < 13.6, 
there exists a finite number of trimer states, the so-called universal trimers (also known as 
Kartavtsev-Malykh trimers l isll ) that carry one unit of angular momentum and have a binding 


energy oc E\E cx a 2 . These states are precursors of the Efimov states occurring at larger 


k2) 
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mass ratios 111 81, 1241. Although they do not exhibit discrete scaling invariance, their existence 
originates from a continuation of the Efimov three-body attraction to mass ratios smaller than 
13.6. For 8.2 < m^/m± < 12.9, there is only one universal trimer state, while there exist two 
universal trimer states for 12.9 < m^/m^ < 13.6. These states are triply degenerate, due to 
their angular momentum £ —1. Unlike the Efimov trimers, the universal trimers are internally 


stable against decay by three-body recombination to possible deeper states [18, 


25] because the 


centrifugal potential between identical']'-']' fermions is strong enough to prevent three particles 
from coming within distances of the order of r 0 where the three-body recombination processes 
can occur. 


According to Ref. [26], a universal tetramer state, composed of three heavy and one light 
particles, exists for m^/m^ > 9.5. Its existence is also related to the occurrence of a four- 
body Efimov attraction for mass ratio m^/m^ > 13.4 [27[]. The existence of larger clusters, 
pentamers or hexamers, is still an open question, but since the Pauli exclusion principle between 
the identical fermions is more significant for larger number of particles, it is reasonable to 
assume that they do not exist in the mass ratio range 8.2 < m^/m^ < 9.5. Therefore, for 
8.2 < m|/?«! < 9.5, the lowest energy few-body bound state is a universal trimer. In this paper, 
we focus on mass ratios around this range, since our main aim is to investigate the possibility 
of a stable trimer phase connected to Efimov physics. 

In order to have a glimpse of the region of parameters in which a gas of trimers is the ground 
state, let us consider the system for a general density imbalance in the limit na 3 — > 0 + limit 
where interactions among trimers, dimers, or atoms can be made arbitrarily weak and the ground 
state is simply the configuration of very weakly-interacting trimers, dimers and unpaired single 
particles. Let us first neglect these weak interactions, which will be studied in the following 
section for the trimer-trimer case, and investigate the phase diagram as a function of density 
imbalance. The key for this analysis is the fact that for the mass ratio 8.2 < m^/m± < 9.5, 
the binding energy of the trimers is larger than, but less than twice, that of the dimer, i.e. 


E, 


( 2 ) 


< 


E, 


(3) 


< 2 


Ei 


( 2 ) 


1181 


2811 . Therefore, the trimers can only appear when the number 


of t fermions exceeds that of | fermions: for n^/n to t < 1/2, the ground state is a Bose- 
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Fermi mixture of dimers and J, fermions as shown in Fig. 2. In the region n^/n tot > 1/2, 
the universal trimers appear, behaving as fermions since they consist of three fermions. For 
1/2 < n-f/n t ot < 2/3, where n tot = + n±, the ground state is a Bose-Fermi mixture of 
dimers and trimers. Since the trimer state carries one unit of angular momentum, the Fermi 
surface becomes triply degenerate, and can be labelled by the rotational number m = ±1,0. For 
a larger population imbalance n^/n tot >2/3, the ground state is a mixture of heavy fermions 
and trimers. 

At n^/ritot = 2/3, there are no dimers or remaining fermions. The ground state in the dilute 
limit, that neglects any trimer-trimer interaction, is a three-component Fermi system of non¬ 
interacting universal trimers. It is remarkable that a three-component Fermi system emerges 
from the original two-species mass-imbalanced Fermi system. From now on we focus exclu¬ 
sively on this region n^/n tot ~ 2/3. 


III. LOW-ENERGY EFFECTIVE HAMILTONIAN OF THE TRIMER PHASE 


Since the centrifugal repulsion between the ^ fermions is strong enough to overcome the 
Efimov attraction at short distances and suppress the 3-body and 4-body Efimov effects in the 
mass ratio window we are interested in 8.2 < m^jm^ < 9.5 L2il23[ 123,1290, we can reasonably 


assume that the same argument applies to 5-body and 6-body Efimov effects. The scattering 
properties between the trimers are therefore universally characterised by the s-wave scattering 
length a between the and | fermions. Each trimer carries one unit of angular momentum 
therefore three scattering channels, labelled by the total internal angular moment F = 0,1, 2 
of the six-body problem, are available. Since the trimers are fermions, the relative orbital state 
has an even parity for F — 1, while it has an odd parity for F = 0, 2. The s-wave scattering 
occurs only in the F — 1 channel, while the dominant channel is the p -wave scattering in the 
F — 0, 2 channels. In the zero-temperature limit, s-wave scattering dominates and the effective 
trimer-trimer interaction is controlled by the F = 1 channel. That implies that interactions 
are restricted to trimers with different projections of the angular momentum. It is important 
to note that the contribution from the other channels is not necessarily negligible in the range 
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of temperatures relevant to experiments. However in line with the exploratory character of 
this study we focus on the zero temperature limit that is well approximated by the following 
trimer-trimer interacting Hamiltonian, 


Hint ~ Y! ^=l( r i - r j) P F=l, (1) 

ij 

where we are assuming that trimers are point particles at r, and r :) , Vp is the corresponding 

F 

interaction potential, P l J = ^ | F, M), 3 , :I (F, M\ is the projection operator onto a two-body 

M=—F 

state of the i-th and j-th trimers in a basis given by the total internal angular momentum F and 
its total projection M. We note that in realistic experimental situations the ground state may not 
only be composed of trimers though this is not important because it is experimentally possible 
to remove other particles from the system. 

Using the Clebsch-Gordan coefficients for F — 1, Eq.([Q) is rewritten in the second quanti¬ 
zation form as 


rr( e ff) _ 

/1 int — 


d 3 r / d 3 r'V F= Ar-r') 


W( r )V\S(OV , o(>')VhM + V’l(r)^- 1 (OV , -i(OVh(r) 


+V’_i(r)^o(r')V’o(OV ; -i(r) 




( 2 ) 


mi m2 


where = ±1, 0 denotes the trimer quantum number related to the projection of the angu¬ 
lar momentum. In the last line above we assume that V F=1 (r — r'), to be computed explicitly 
in the next section, is a short-range potential. This assumption makes the terms correspond¬ 
ing to mi = m 2 negligible due to Fermi statistics. Interestingly, the last line in Eq. © is 
SU(3) symmetric [30], even though the initial heavy-light mixture has much less symmetry. 
This enhanced SU(3) symmetry could be observed experimentally provided that the potential 
is sufficiently short-ranged and other processes that break the degeneracy of m-states, such as 
dipole-dipole interactions, are negligible. 

Finally we note that due to the short-range form of Eq. ©, spin-flipping processes are 
heavily suppressed. Therefore, the initial population in each internal angular momentum state 
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n±i,n 0 is an almost conserved quantity. In order for the system to evolve into an equilib¬ 
rium state n_i = ri\ = n 0 , it is necessary to consider spin-flipping processes occurring in the 
higher-partial-wave scattering channels. We shall see later that this sets a natural time scale for 
equilibration. 


IV. TRIMER-TRIMER POTENTIAL AND THE NATURE OF THE GROUND STATE 


The nature and stability of the trimer phase strongly depends on the trimer-trimer interaction 
potential V F=1 . If it is purely repulsive, then > 0 and the ground state is a three-component 
Fermi liquid composed of trimers. If it is sufficiently attractive so that a t I t : =1 < 0, then it leads to 
an SU(3) superfluid state [13 lL 1321 in which the building blocks are not pairs of particles but pairs 
of trimers. However, if the attraction is too strong, this trimer phase can be unstable against col- 
lisional losses in ultracold atom experiments, which occur when the atoms of colliding trimers 
come within distances of the order of r 0 and recombine into more deeply bound molecules. 

To understand better the trimer phase, we need to carry out an explicit calculation of the 


trimer-trimer potential and the associated s-wave scattering lengt 


is reminiscent of the problem of the baryon-baryon potential [33 


l in the F = 1 channel. This 


35] in nuclear physics within 


the quark model. In both problems, basic constituents are fermions, so the Pauli principle is 
expected to play an important role especially for short distances. However, there are important 
differences. For instance, in the universal trimer case, a long-range attraction between the heavy 
t fermions can be induced by exchanging the light,[ fermions between the trimers. This is in 
marked contrast to the short-range character of typical nuclear forces. The determination of 
the trimer-trimer interaction is a rather non-trivial problem, requiring a six-body calculation. 
Solving exactly the six-body problem, however, is a very demanding numerical calculation, so 
we opt here for approximate methods. 

We first compute the trimer-trimer potential by the Bom-Oppenheimer (BO) approxima¬ 
tion. The three-body physics, including the appearance of the universal trimers and the Efimov 
trimers, can be adequately reproduced by the BO approximation thanks to the large mass imbal¬ 
ance. It is therefore natural to extend the BO approximation to the six-body problem as a simple 
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but heuristic approach to understand the trimer-trimer scattering process. In the BO approxi¬ 
mation, the light 4 , fermions’ Schrodinger equation is solved by fixing the heavy f particles’ 
positions. The eigenenergy of the light fermions acts as an adiabatic potential between the t" 
particles. Subtracting the internal energy of each trimer, assuming the two pairs of heavy par¬ 
ticles are well separated, we obtain the trimer-trimer potential (see Supplementary information 
for more details [50]). In Figs. 3 (a) and (b), we show the trimer-trimer potential calculated with 
the BO method for various heavy particles’ configurations sketched in Figs. 3 (c) to (h). The 
trimer-trimer potential is repulsive at large distance for all configurations, decaying exponen¬ 
tially as the dimer wave function of the light fermion around the heavy particle. The repulsion 
is due in part to the fermionic nature of the light particles: although the first light fermion oc¬ 
cupies a bonding orbital, the second light fermion has to occupy an antibonding orbital which 
overcomes the attraction of the bonding orbital. We note that the centrifugal repulsion between 
the heavy fermions are not included in the BO potential shown in Fig. 3. Therefore, the actual 
trimer-trimer interaction would be more repulsive than that shown in Fig. 3. 

The BO approach has a several shortcomings: it is not realistic to compute the trimer-trimer 
scattering length as it requires solving a 4-body Schrodinger equation with the long-range BO 
potential, which would be hardly easier than solving the original six-body problem with a con¬ 
tact interaction. Moreover, it cannot correctly capture the trimer-trimer potential at short dis¬ 
tance where, see Fig. 3, the potential is attractive for most of the configurations. We note that 
the attraction is likely an artefact of the way in which we subtract the trimer internal energy, 
which is strictly valid only when the two trimers are well separated (see Supplementary infor¬ 
mation [50]), and the fact that the centrifugal repulsion among heavy fermions is not included. 
In summary, the BO results depicted in Fig. 3 could be interpreted as the trimer-trimer potential 
only when the trimer-timer separation is much larger than the trimer size R> a. 


Faced with these limitations, we use a second approximate method to evaluate the trimer- 
trimer potential and the scattering length that addresses satisfactorily some of these problems: 
the resonating group method (RGM) [36-39] commonly employed in nuclear physics calcula¬ 
tions. Unlike the BO method, the RGM takes into account exactly the antisymmetrization of 
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the wavefunctions due to the Pauli principle which typically plays an important role for short 
distances. The RGM approximates the six-body wave function describing the trimer-trimer 
scattering by the antisymmetrized product of the known trimer wave functions and an unknown 
wave function describing the relative motion between the two trimers. Therefore, it is a good 
approximation provided that, except for the exchange of identical particles, the wave functions 
for the three-body subsystems are not substantially altered during the scattering process. More 
specifically, the RGM approach is qualitatively correct unless the trimer-trimer interaction is 
strongly attractive so that there is a six-body bound state. Although we cannot rule out that this 
is the case, the fact that a tetramer state only appears at a larger mass ratio [26] than the ground 
state trimer suggests that, a six-body bound state, if exists, it must occur at an even larger mass 
ratio. 


From the RGM ansatz, one can derive an effective interaction potential between the two 
trimers [40]. As is shown in Fig. 4 (a), we have found that this potential is repulsive for all 
distances. For large separations, the decay of the potential is exponential, in agreement with 
the BO results. In Fig. 4 (b), we present the trimer-trimer scattering length for different mass 
imbalances obtained from the RGM trimer-trimer potentials. For all mass ratios between 8.2 
and 12.9, for which only one trimer state exists, the scattering length is found to be positive. 

Although not exact, the RGM method indicates that, due to Fermi statistics, trimers expe¬ 
rience a repulsion for separations of the order of the scattering length and larger. The RGM 
prediction that the interaction is still repulsive at smaller separations is valid provided that the 
trimer-trimer scattering does not result in a six-body bound state in the range of mass ratios 


of interest. This strongly suggests the possibility that colliding trimers, like dimers [41, 


m, 


are protected from collisional losses occurring at shorter distances, raising the prospect of the 
experimental realization of the trimer phase in ultracold-atom experiments. As to the nature of 
the trimer phase, the positive trimer-trimer scattering length obtained in our RGM calculation 
indicates that it should be an SU(3) Fermi liquid, rather than an SU(3) superfluid. The RGM 
calculation being approximate, we cannot completely rule out the possibility of the latter. In¬ 
deed, we expect that an important correction, that increases with the mass ratio, comes from the 
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virtual excitations of the trimers into the dimer-particle continuum. Perturbation theory to sec¬ 
ond order indicates that this correction is attractive. An explicit calculation would be required 
to clarify whether, for certain mass ratios, the strength of this correction is strong enough to 
destabilize the Fermi liquid ground state. The most interesting scenario would correspond to a 
small attraction at large distance inducing superfluidity, while preserving the stabilizing Fermi 
repulsion at short distance. 


V. EXPERIMENTAL REALIZATION 


Prime candidates to observe experimentally this novel trimer phase are fermion mixtures 
of ultracold atoms. A mixture of 6 Li and 53 Cr atoms [43] (rn^/m^ = 8.80) falls within the 
mass-ratio window 8.2 < m^/m^ < 9.5 [44] in which a trimer Fermi liquid phase is possible. 


If the s-wave scattering length is tuned to satisfy the condition A, n 


-1/3 

tot 


» a > r vdw via a 


Feshbach resonance, where r vdw is the van der Waals length of the atoms, the trimer Fermi 
liquid phase can be realized at the density imbalance /n tot & 2/3. Another experimental 
setting of interest is a Li-K mixture where radio-frequency techniques have shown[45] a strong 
atom-dimer p-wave attraction at a lower mass ratio of 6.63. 

In any case a necessary condition for the observation of the trimer phase is that the typical 
time scale for its formation and equilibration is within the reach for current cold-atom experi¬ 
ments. Without any external influence, the trimer phase forms via a four-body recombination 
of the fermionic mixture, which is too slow to be observed in cold-atom experiments. However, 
trimers can be formed within a much shorter time compatible with experimental observation 


by using an adiabatic change of the s-wave scattering lengf 


frequency association of trimers in an atom-dimer mixture [48 


l 147], or performing the radio- 
49]. Once the trimers are formed 


by these methods, one must wait for the trimer gas to equilibrate. There are two important time 
scales for the equilibration: thermalization time r th and relaxation time of the trimer’s internal 
states Tint- Since the trimers can change their internal states only via the p-wave collision, the 


relaxation time r re i ~ h(k 


( tot ) a )- 6 


3 /4 tot) , where kf l) 


and E 14 are the total Fermi wavevec- 


tor and Fermi energy of the mixture respectively, is much larger than the thermalisation time 
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r t h ~ h(k ( F 0t) a)- 2 /E^ > , which occurs via s-wave scattering. Therefore, the relaxation time 
gives the minimum required time in experiments to realize an equilibrium trimer Fermi liquid 

(tot) -■ O.OOSnm^ 1 , corresponding to a density n ~ 10 18 atom/m 3 , 

0.3, r t h ~ 0.1 ms 

and r re i ~ 10 ms, which is an accessible time scale in cold-atom experiments. 

Another potential source of deviations between theory and experiments is the presence of 
magnetic dipole interactions in the mixture that break the 577(3) symmetry of the model. This 
effect will be noticeable only if the energy related to the magnetic dipole interaction E D ~ 
/r 2 /io/47rr 3 is comparable with the trimer bound energy Eb ~ ti 2 /2m 2 a 2 . Assuming r ~ a, 
as the typical size of the trimer, /i 0 the vacuum permeability and // « 6/j r with the Bohr’s 


-,(tot) 


phase. For a r^i lOOnm , k F 
Ep 0t) = h 2 kp 0t ^ /2m, with m here the mass of the atom, we obtain kp 0t ^a ~ 


magneton [46] it turns out that, for any a lnm, a condition easily met experimentally, the 


splitting energy Eb is always much smaller than Ed- This is a strong indication that the ,57/(3) 
symmetry is to be approximately conserved in realistic experimental conditions. 


VI. CONCLUSION 

We have identified a region of parameters in a mass-and-density-imbalanced two-species 
fermionic mixture for which the system undergoes a transition to a weakly-repulsive SU(3 ) 
Fermi liquid composed of stable universal trimers. The parameters required to realize this 
many-body phase are within reach with cold-atom experiments on mass-imbalanced fermionic 
mixtures, such as 6 Li and 53 Cr atoms. 

The present work also raises many interesting questions about the role of Efimov physics 
in many-body systems. What happens in the resonantly-interacting regime n to ta 3 —* oo ? 
In which circumstances can Efimov physics induce superfluidity? Is the ground state for larger 
mass ratio rn^jrn^ > 9.5 a Bose gas of universal tetramers [26] ? How do the trimer and tetramer 
phases evolve with mass ratio and population imbalance? The answers to these questions would 
be an important advance in the understanding of the conditions for Efimov-related clusters to 
constitute the building blocks of novel phases of quantum matter. 
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FIG. 1. The number of f (n-j-) and 4 , (njJ fermions, dimers n ( \\m and trimers n tr j as a function of 
the population imbalance n^/ntot, where ntot = n -\ + n i is the total fermion density. These results are 
strictly valid in the dilute limit that neglects completely interactions among fermions, dimers and trimers. 
Moreover for fermionic ground states we also neglect the effect of a finite Fermi pressure. In this paper, 
we focus in the region around n|/n l0 i = 2/3 where the ground state is a superposition of universal 


trimers. 
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FIG. 2. (a) (b) Trimer-trimer potential, in units of h fir no 1 with rn the mass of the light particle, calcu¬ 
lated with the Born-Oppenheimer approximation for various configurations illustrated in (c)-(h). From 
(c) top (in the large R/a limit) red curve to (h) bottom cyan (light grey in a greyscale) curve. The dis¬ 
tance between the heavy particles is taken to be q = 1.6a. The trimer-trimer interaction is repulsive at 
large distance for all configurations. It decays exponentially with a decay constant very similar to that 
of the square of the dimer wave function (black dashed line in Fig. (b)) corresponding to a light particle 
around a heavy particle. 
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(a) 



(b) 



FIG. 3. (a) Trimer-trimer potential for rrq/rn; = 9 calculated with the resonating group method, where 
//f; = m|/2 + m| is the reduced mass between the trimers. It is repulsive for any trimer-trimer separation 
R and azimuthal angle 6 between the trimers measured from the quantization axis of the trimer internal 
states, (b) Trimer-trimer scattering length a^ =1 as a function of mass ratio calculated with the RGM 
method. It is always positive and decreases as the mass ratio increases. This is expected as the trimers 
are more tightly bound for larger mass ratio and therefore then - cross section is smaller. We note that, 
strictly speaking, our results are only applicable for mass ratios < 12.9. For larger mass ratios an 
additional universal trimer state becomes available and the calculation of the scattering length becomes 


much more involved. 




